Solid-State Lasers for Coherent Communication 
and Remote Sensing 

NASA Grant NAGW-1760 


Semi-Annual Progress Report 
for the period 

October 1,1989 through March 31, 1990 

Robert L. Byer 
Principal Investigator 


Edward L. Ginzton Laboratory 
W.W. Hansen Laboratories of Physics 
Stanford University 
Stanford, California 94305 


June, 1990 



Solid-State Lasers for Coherent Communication and Remote Sensing 


NASA Grant NAGW-1760 


Robert L. Byer 
Applied Physics Department 
Stanford University 
Stanford, California 94305 


ABSTRACT 


Laser development, high efficiency, high power second harmonic 
generation, operation of optical parametric oscillators for wavelength diversity and 
tunabihty , and studies in coherent communications are reviewed as achievements of 
the second semi-annual performance period of this program. 
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I Introduction 

NASA contract Grant NAGW-1760 supports development of solid state lasers and 
nonlinear frequency conversion of these sources for remote sensing and coherent optical 
communication applications. In the following pages, progress made during the last half 
year( October 1 1989 to March 31 1990) is reported. Direction for future research is also 
indicated. 

II Review of Progress 
A. Harmonic Conversion 

During the last half year, emphasis has been placed on extending previous work on 
efficient single frequency externally resonant Second Harmonic GenerationfSHG)! 1 ! to the 
high power regime. Watt level coherent single frequency green radiation is ideal as a laser 
source for gravitational wave observatories. Other applications of the high power SHG 
source include the use as a pump source for coherent nonlinear optics studies, and when 
the 532 nm radiation is further up converted into the UV, it can be used as a light source for 
fine line VLSI lithography. 

Previously, theoretical calculations have shown that multi-watt 532 nm radiation 
can be produced in an external cavity configuration with a 10 watt 1.06 |im pump using a 
number of different nonlinear optical crystals. These crystals include: BaB^CU, LiB 3 <D 5 , 

[1] W.J.Kozlovsky,C.D.Nabors and R.L.Byer,"Efficient Second Harmonic Generation of a Diode-Laser- 
Pumped cw Nd: YAG Laser Using Monolithic MgO:LiNbC>3 External Resonant Cavities," IEEE J. 

Quantum Electron. 24, pp.913-919(June,1988). 
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and KTP. The more traditional Mg 0 :LiNb 03 crystal is not included due to its relatively low 
photorefractive damage threshold. Among the crystals examined, UB 3 O 5 appears most 
promising since it can be 90 degree phase matched at 1.06 pm wavelength, and it is 
reported to be very tolerant of high laser intensity. 

In order to facilitate evaluation and comparison of the different crystals, we have 
opted for a discrete rather than monolithic resonator during initial trials. The enhancement 
resonator we have constructed has a bow tie configuration. It consists of two flat mirrors 
and two curved mirrors. The nonlinear crystal is placed between the two curved mirrors 
where a tight focal spot is located. One of the flat mirror is mounted on a piezo electric 
mount so that the cavity resonant frequency can be servo locked to the laser frequency. The 
pump source for the high power harmonic experiment is the injection locked 13 Watt 
Antares laser that has been previously developed in the group. f 2 l 

So far, we have tried both KTP and BBO crystals. 532 nm radiation of 500 mw 
and slightly more than 1 watt has been produced using the two crystals respectively. These 
are less than what we expected, and the limiting factor right now is the excessive single 
pass power absorption of the crystals. We are currently setting up a very sensitive 
calorimeter for carrying out detailed crystal absorption measurements. By measuring the 
temperature rise in the center of the crystal, we can also tell if the absorption is due mainly 
to surface scattering or bulk absorption. 

In addition to characterization of crystal absorption, we are also continuing the 
measurement of nonlinear optical coefficients of different crystals using the phase matched 
second harmonic technique. These measurements were carried out initially using a pulsed 
laser. The results of these measurements have been written up for publication. A draft of 
the paper is included in appendix A. At the time when the initial measurements were made, 
several discrepancies with previously published values were found. In addition, 

[21 C.D.Nabors, A.D.Farinas, T.Day.S.T.Yang, E.K.Gustafson, and R.L.Byer/’Injection locking of a 13- 
watt NdrYAG ring laser," Opt. Lett. 14, 1189-1191(1989). 
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measurements of nonlinear coefficients of LiIC>3 using SHG technique and parametric 
fluorescence technique disagree.To corroborate the pulsed measurements, CW SHG 
measurements were subsequently carried out using the injection locked Antares laser. In 
general, the results from the CW measurements have agreed with the pulsed results. 
However, disagreements with parametric fluorescence measurements of UIO3 remain a 
puzzle. We are continuing to investigate the cause for such discrepancies. 

B. Optical Parametric Oscillator(OPO) 

1 . Doubly Resonant OPO (DRO) 

Since the demonstration of the single axial mode operation of CW Doubly Resonant 

OPO(DRO) in our group, t 31 we have concentrated on investigating the coherence and 

tuning properties of CW DRO. The most significant result from the coherence studies is the 

confirmation that the DRO output preserves the pump coherence. The linewidth of the 

beatnote between the DRO signal mode and an independent single frequency Nd:GGG 

laser shows the same 13 KHz width as the beatnote between the Nd:YAG pump and the 

Nd:GGG laser, indicating that the DRO adds negligible noise to the pump. Furthermore, 

DRO output is found to phase locked to the pump when operated at degeneracy. Thus, at 

degeneracy, the DRO produces a phase coherent subharmonic of the pump. This has 

implication for application of the DRO as an optical frequency divider. Away from 

degeneracy, the phase of the signal and idler are found to be anti-correlated with each other. 

The results of these experiments have been published. A reprint is included as appendix B. 

DRO tuning and mode control has also been studied during the last half year. Due 

to the requirements for simultaneous oscillation of signal and idler mode and the 

conservation of energy in a DRO, DRO operation is very sensitive to pump frequency 

[3] C.D.Nabors, R.C.Eckardt, W.J.Kozlovsky, and R.L.Byer, "Efficient, single-axial-mode operation of a 
monolithic Mg0:LiNb03 optical parametric oscillator," Opti. Lett., 14 , 66-68(1989). 
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fluctuation, cavity length changes and crystal oven temperature drift. By using a stable 
single frequency diode pumped solid state laser as the pump and resorting to the monolithic 
cavity design, the effects of the first two perturbation are reduced. It is now possible to 
contemplate stabilized single frequency operation and continuous tuning of the DRO by 
changing the crystal temperature to adjust phase matching requirement, vary pump 
frequency to adjust energy conservation relationship, and then to tune the cavity length by 
applying electric field to the crystal and thus modified the simultaneous resonant conditions 
for signal and idler modes. To understand the control necessary to achieve slow tuning and 
stabilized single frequency operation of DRO, an extensive theoretical analysis has been 
earned out using the Mg 0 :LiNb 03 monolithic DRO as a specific example. By taking into 
account the temperature dependent crystal dispersion, crystal thermal expansion, the electro 
optic effect and the piezoelectric effect, the theory can be used to predict the tolerance 
required on pump frequency stability, temperature drift, and cavity length change to achieve 
stabilized single frequency operation and slow tuning of DRO. Initial experimental results 
using the Mg 0 :LiNb 03 monolithic DRO show close agreement with theory. Additional 
experiments are under way to demonstrate tuning and stabilized single frequency operation 
of DRO. 

For these experiments, a novel 3: 1 DRO made of congruent LiNb 03 will be used. 
When pumped with 532 nm radiation, the 3:1 DRO will produce signal and idler 
wavelengths around 2 micron and 700 nm wavelengths respectively. The 2 micron regime 
is particularly interesting since it contains several water absorption lines. The 3:1 DRO can 
thus be used for remote monitoring of moisture in air. In addition, 2 micron radiation is at 
an eye safe wavelength and it is interesting for use for wind velocity LIDAR. Beside 
remote sensing and LIDAR application, the 3:1 DRO will serve as a 3 to 1 frequency 
divider that will be useful as part of the frequency divider chain that links the Cesium 
atomic clock reference to other parts of the spectrum. 
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2. Singly Resonant OPO (SRO) 


Compared to DRO operation, SRO operation is not as highly constrained by requirement 
for simultaneous resonance of signal and idler. As a result, continuous tuning is much 
simpler in SRO than in DRO. For this and other reason, demonstration of CW SRO 
operation has been a long standing goal in our group. Unfortunately, previous attempts to 
achieve CW SRO operation were unsuccessful due to material damage limitations. Since 
then, we have started investigation of other suitable nonlinear optical materials for use as 
SRO. As our next attempt, we will investigate the transition from DRO to SRO operation 
by tuning the 3:1 DRO operating frequency off the mirror coating range. The pump source 
will be either a single line 5 Watt Ar ion laser or the 4-5 watts produced from the external 
resonant doubler. If the initial trial is successful, application of SRO as a tunable chirped 
optical radar will then be studied. 

3. Squeezed State Generation 

Due to its unique quantum mechanical properties, OPO has been considered a prime 
candidate for squeezed state generation. As pan of our OPO research, we are also interested 
in investigating its use as a practical squeezed state generator. Applications of such 
squeezed state generator include ultra precise interferometric measurements such as what is 
required for gravitational wave detection. In addition, application of squeezed state for 
optical communication is also been actively pursued. 

One possible way to generate “squeezed” state radiation from OPO is by taking 
advantage of the close correlation between signal and idler photons resulting from the 
parametric process. Heidmann et.al. in France has earned out an elegant experiment that 
demonstrates application of signal idler correlation for sub-shot noise detection.! 4 ! In this 

[4] A.Heidmann, RJ.Horowicz, S.Reynaud.E.Giacobino, and C.Fabre,"Observation of Quantum Noise 
Reduction on Twin Laser Beams," Phys. Rev. Lett., vol. 59,no.22, pg. 2555-2557,1987. 
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experiment, the signal and idler intensities are separately detected on two nominally 
identical detectors. Due to signal and idler correlation, the subtracted photo current from the 
two detectors falls below the vacuum shot noise level within the bandwidth of the OPO 
cavity. In their work. Type II phase matching was used which allows for separation of 
signal and idler modes by polarization discrimination. 

In our work, we attempted to demonstrate the sub-shot noise detection by using 
type I phase matching in a monolithic OPO. To separate the signal and idler modes, the 
OPO is forced to operate very far away from degeneracy so that the wavelength difference 
between the two is large enough to allow separation by prisms. The best result obtained 
shows about 50% noise reduction of the subtracted photo currents compared to the vacuum 
noise level. The noise reduction is presently limited by detector efficiency, and absorption 
loss in the OPO crystal. The results of the squeezed state experiment has been submitted for 
publication. A pre-print is included in appendix C. 

As a demonstration of an application for such a “Two Color squeezed” light, we 
injected the signal and idler beams into the two arms of a Mach Zehnder interferometer, in 
which one of the arm has a Pockel cell that imposes a slight amplitude modulation. We 
show in the experiment that sub shot noise amplitude modulation can be observed by 
looking at the subtracted photo currents between signal and idler. This has direct 
application for laser spectroscopy in detecting weak absorption lines. 

C. Application of Diode Pumped Solid State Laser in Coherent 
Communication 

The advent of diode pumped solid state laser dramatically reduced the inherent linewidth of 
solid state lasers. Presently, beatnote linewidth of 10 to 20 KHz can be routinely obtained. 

By locking to an external reference via electronic servo technique, the linewidth of these 
lasers can be further reduced. A beatnote linewidth of 300 mHz between two NonPlanar 
Ring Oscillators(NPRO) locked to a high finesse optical cavity has recently been 
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demonstrated in our laboratory.^ This work was sponsored by NASA SUNLITE (NAG1- 
839) program. Applications of such a highly stabilized laser include gravitational wave 
detection, high resolution laser spectroscopy, and coherent communication. 

Under the present program, we are particularly interested in exploring the 
application of low noise NPROs for coherent communication applications. As a laboratory 
demonstration of a communication link using Phase Shift Keying(PSK), an Optical Phase 
Lock Loop(OPLL) has been set up using two NPRO lasers serving as the transmitter and 
local oscillator. The details of this work are included in appendix D. By locking the phase 
of these two lasers at quadrature, any phase modulation imposed on the transmitter can be 
demodulated by looking at the error signal of OPLL. Coherent communication using 
homodyne PSK is the most sensitive data transmission technique. Semiconductor lasers 
have not been able to take advantage of extra sensitivity offered by homodyne PSK due to 
their inherent large phase uncertainty. Previous attempts to realize OPLL using 
semiconductor lasers have involved elaborate schemes for narrowing the linewidth of these 
lasers. Due to inherent low noise, diode pumped solid state lasers are ideal for OPLL 
applications. Indeed, a modest bandwidth servo loop proved sufficient to phase lock the 
two NPRO lasers with less than 12 mrad rms phase error. Theoretical analysis of minimal 
phase error in an OPLL also brings out additional advantages of diode pumped solid state 
lasers over semiconductor lasers in homodyne communication links. It is found that since 
the linewidth of a diode laser is dominated by broad band quantum noise contributions 
whereas the diode pumped solid state laser linewidth is low frequency flicker noise 
dominated, residual phase noise error of OPLL can be lower for diode pumped solid state 
lasers versus diode lasers with the same linewidth. This conclusion also implies that low 
data rate communication is feasible using diode pumped solid state lasers without incurring 
excessive crosstalk penalty. In our laboratory demonstration system, data rates as low as 
20 KHz can be used. This is an important consideration especially for a space based 
f 5 l T.Day, Phd. dissertation, Stanford University, 1990. 
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coherent communication link where data rate and laser power requirements can be traded 
off against each other. 

Ill Future Research 


In the next half year, investigation into coherent solid state laser sources for remote sensing 
and optical communication applications will continue as outlined in the original proposal. In 
the next quarter, emphasis will be placed on improving the conversion efficiency of discrete 
cavity external SHG. This will involve better characterization of crystal absorption using a 
calorimeter that is capable of 0.01% absorption measurement. In addition, investigation of 
novel nonlinear optics materials will continue as they become available. 

Research in OPO will concentrate on demonstration of stabilized single frequency 
operation and slow controlled tuning of the 3:1 OPO. Transition into the SRO regime will 
also be investigated using the same 3:1 OPO by tuning off the mirror reflectivity range. 
Longer term research goals include demonstration of chirped optical r adar using an SRO 
and application of the DRO for an optical frequency divider chain. Finally, more 
speculative research into squeeze state generation using OPOs will continue. 

In the area of solid state laser development for coherent communication, the 
research direction will be focused on higher power, low noise, single frequency, diode 
pumped slab laser to serve as local oscillator in a space based coherent communication link. 
The laser under development for gravitational wave astronomy will satisfy all the 
requirements for optical communication. The proposed slab laser will be pumped by sixty 1 
watt, fiber- coupled single stripe diode laser to produce 10 watts of output. The slab laser 
will then be injection locked to achieve single frequency operation, and stabilized to a high 
finesse optical cavity to reduce the linewidth further. With the progress that has been made 
and the planned experiments yet to be carried out, we fully expect to significantly improve 
solid state laser technology and nonlinear frequency conversion of these lasers for remote 
sensing and optical communication applications. 
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Absolute and relative nonlinear optical coefficients of KDP, KD*P, BaB 2 0 4 , LiI0 3 , 
MgO:LtNb0 3 , and KTP measured by phase-matched second harmonic generation 


Robert C. Eckardt, Hisashi Masuda, Yuan Xuan Fan, and Robert L. Byer 


abstract 

Both absolute and relative measurements of the nonlinear optical coefficients of 
six nonlinear materials were measured by second harmonic generation. A 
single-mode, injection-seeded, Q-switched Nd:YAG laser with spatially filtered 
output was used to generate the 1 .064 pm fundamental radiation. The following 
results were obtained: r/ 3 g(KDP) = 0.38 pm/V, P 3 ^(KD*P) = 0.37 pm/V , 
W 22 (BaB 2 0 4 )l = 2.2 pm/V, J 31 (LiI0 3 ) = -4.1 pm/V, r/ 31 (5%MgO:LiNbC> 3 ) = 
-4.7 pm/V, and d eff <KTP) = 3.2 pm/V. The accuracy of these measurements is 
estimated to be better than 10%. 


This work was supported by U. S. Army Research Office contract DAAL03-88-K-01 13 and NASA Grant 
NAGW 1760. 

R. Eckardt and R. Byer are with the Ginzton Laboratory, Stanford University, Stanford CA 94305. 

H. Masuda was a visiting scicnust at the Gin/ton Laboratory, Stanford University, Staniord CA 94305. He 
is with the Sony Corporation, Corporate Research Laboratories, 1-7-4, Kohnan, Minato-ku, Tokyo 108, Japan. 

Y. X. Fan is with LaserGenics, Inc., San Jose, CA 95131. 
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Absolute and relative nonlinear optical coefficients of KDP, KD*P, BaB 9 0 4 , UI0 3 , 
MgO:LiNb0 3 , and KTP measured by phase-matched second harmonic generation 


Robert C. Eckardt, Hisashi Masuda, Yuan Xuan Fan, and Roben L. Byer 


I. Introduction 

Careful measurements of the nonlinear optical coefficients of several nonlinear 
crystals of current technical interest are presented. These measurements were made by 
phase-matched second harmonic generation using a Q-switched Nd.YAG laser operating at 
1.064 (im as the source of fundamental radiation. The laser was injection seeded to 
operate in a single longitudinal mode, and the output was spatially filtered to produce a 
Gaussian-like transverse distribution. Characterization of the pump pulse and calibration 
of the apparatus indicated that these measurements should be accurate to better than 10%. 

This measurement was originally intended to be only a determination of the nonlinear 
coefficient of barium metaborate (BaB 2 0 4 ), an interest that was motivated by observations 
of optical parametric oscillator thresholds that were below expected values [ 1 ]. The 
measurement grew to include several other nonlinear optical materials that are commonly 
used for second harmonic generation pumped by 1.06 pm radiation. The additional 
materials were measured because of a current controversy concerning the absolute 
magnitude of nonlinear coefficients. 

Our observations use both absolute and relative measurements of nonlinear 
coefficients. The absolute measurements yielded values consistent with other second 
harmonic measurements. Measurements performed with high-power well characterized 
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lasers and potassium dihydrogen phosphate [2] have yielded the value c/^(KDP) = 0.39 
pm/V, which has become accepted in the application of high-power harmonic conversion 
for fusion research [3]. Earlier cw measurements of 633- to 316-nm harmonic generation 
in ammonium dihydrogen phosphate gave a value of = 0.57 pm/V [4], [5]. 

Combining this with relative measurements between KDP and ADP yields t/ 36 (KDP) = 
0.41 pm/V [6]. We measured t/ 36 (KDP) = 0.38 pm/V in good agreement with both earlier 
harmonic measurements. Values for </ 36 (KDP) determined by the technique of parametric 
fluorescence [7], [8] tend to be higher. Previously the value for 1.06 pm to 532 nm 
harmonic generation in lithium iodate t/ 3 j(LiI0 3 ) = -7.1 pm/V based on parametric 
fluorescence has been accepted as standard [9]. (No attempt was made to measure signs 
of the nonlinear coefficients, however when the signs are known [6] they are included.) 

On the parametric fluorescence LiI0 3 scale, the coefficient tor potassium dihydrogen 
phosphate is d 36 (KDP) = 0.63 pm/V for second harmonic generation pumped by 1.06 pm 
radiation [7]. Clearly there is a discrepancy in the accepted values for nonlinear optical 
coefficients. 

The crystals included in the measurements reported here are KDP, LiI0 3 , BaB 2 0 4 , 
potassium dideuterium phosphate (KD*P), potassium titanyl phosphate (KTP), and 
magnesium oxide doped lithium niobate (5%Mg0:LiNb0 3 ). The ratio ot t/ 36 tKDP) to 
d 31 (LiI0 3 ) obtained in these measurements is in agreement with earlier measurements [7]. 
However, the value of d 31 (LiI0 3 ) is only 58% of the parametric fluorescence value 
presented in [9]. The reason for the variance is not understood. It has been suggested that 
the pump wavelengths may have been too close to the LiI0 3 absorption edge in the 
parametric fluorescence measurements. There has also been some question of the 
necessity for further corrections to deal with the birefringence of the nonlinear crystals. 

An anisotropic Green's function analysis [10] is discussed in the appendix with the 
conclusion that existing theory is appropriate for the quantitative analysis ot near field 
second harmonic generation and parametric fluorescence. 
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The advantage of parametric fluorescence is that it is only necessary to measure the 
power ratio of the pump radiation and the generated fluorescence, whereas second 
harmonic requires the measurement of an absolute power and the spatial an temporal 
distributions of the pump field. Phase matching is also of critical importance in the second 
harmonic technique. Second harmonic generation, however, has the advantage that it is a 
direct measurement of the nonlinear coefficient under conditions that more closely 
duplicate actual applications. The results presented here are a study of nonlinear 
coefficient measurement by second harmonic generation. At this point no attempt is made 
to explain the nonlinear coefficient values obtained by parametric fluorescence compared 
to those obtained by second harmonic generation. However, clearly more work is 
indicated in the future regarding the nonlinear coefficient values determined by the 
fluorescence method. 

Our second harmonic measurements were performed under near field conditions 
where diffraction is not of concern. Birefringent walkoff, however, was significant in the 
longer angle-tuned crystals. A Gaussian transverse intensity distribution was assumed, 
and the focused beam analysis of Boyd and Kleinman [11] was used. This cw analysis 
was extended to pulsed harmonic generation by numerical integration of the observed 
temporal distribution. Substantial effort was placed in characterizing the temporal and 
spatial distribution of the pump pulse. The largest single source of error in these 
measurements is fluctuation in spatial distribution caused by drifting of spatial filter 
alignment This single error source is about equal in magnitude to all other sources of 
error combined. 

The highest levels of harmonic conversion used here were 14%, a level at which 
pump depletion must be considered but can be treated adequately with a simple near field 
approximation. The depletion approximation was verified by observation of the 
dependence of harmonic conversion on pump energy in two crystals. Measurements to 
demonstrate that two photon absorption was not a factor were also performed in the same 
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crystals. In all cases phase matching tuning curves were carefully compared with 
calculated curves obtained from dispersion relations. The agreement of observed and 
theoretical tuning curves assured that the full length of the nonlinear crystal was used, and 
that optical distortions were not limiting harmonic conversion. 

The theoretical development used to analyze our second harmonic measurements is 
discussed in section II. Although no new theory is presented, except for the depletion 
approximation, the theory is presented to make explicit the definitions that are used to 
describe the second harmonic generation interaction. Different factors appear in the 
various methods of presentation, and misinterpretation of factors abound in this field. The 
experimental apparatus is described in section III, and measurements are discussed in 
section IV. Results are summarized in section V. 

II. Theory 

The analysis used here draws almost entirely from earlier work [11], [12] which is 
reviewed in many places [6], [9], [13], ( 14], MKS units are used. Definition of notation 
is made using a few expressions appropriate for monochromatic planewaves. The 
discussion continues by including the focused beam analysis and is extended to pulsed 
harmonic conversion by numerical integration using observed pulse shapes. A simple 
approximation derived empirically by calculation is used to allow tor the small amount of 
fundamental wave depletion in the measurements. Finally a near field approximation is 
used to allow for an elliptical transverse beam distribution. 

The electric field or electric polarization, both of which are real, can be expressed as 
the product of a complex amplitude and exponential summed with the complex conjugate 
of that product. In this format, the time dependent electric field of a monochromic 
planewave of angular frequency to is 
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E ( r,t ) - -i E (r, co ) exp{i (k r - m )) + c.c. 

The relationship between the vector components of electric polarization at the harmonic 
frequency 2 co resulting from the second order nonlinear interaction define the nonlinear 
optical coefficients d ikl (-2co,co,co) 

3 

NL V 

& (r, 2(0) = 2^£ 0 d i}k (-2co,oj,co )£.(r, co )£ k (r, co ) , (1) 

j,k = l 

where % is the permittivity of free space. The reduced notation d i]k — > d m allows the 
representation of the nonlinear optical coefficients in the customary 3x6 matrix. Further 
reduction to a single effective nonlinear coefficient d cff dependent on nonlinear optical 
properties, phase matching, and crystal orientation is standard for modeling phase- 
matched second harmonic generation. With this notation the coupled equations describing 
harmonic generation tor a monochromatic planewave propagating in the z-direction are 

d „ . -itXk.2 „ 2. 

-^E(z,2co) = i kq E (z, co) and (2a) 

d i /\ it z * 

-^E(z,co) = i kz E(z,2co)E (z,co) , (2b) 

where K= co d^f (n c ) with n the index of refraction, c the speed of light, and 
Ak = *2to ' the wavevector mismatch . Two elementary solutions to the coupled 
equations are 

/ 2(i) (0 = IJ.0){n sin{Akl / 2) / (Akl / 2)} 1 and (3a) 

/(/)=/ (0) tanh 2 (H ) (3b) 

ZCO CO 

for the respective cases of negligible pump depletion and perfect phasematching. Here 
/ 2a) (/ ) is the harmonic intensity generated in a crystal of length /, /^(O) is the initial 
fundamental intensity, and T 2 = 2 k 2 !^) / (c n e Q ). Any of Equations (1), (2) or (3) 
can serve to define the notation used to express the nonlinear coefficients. 
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The dispersion in the nonlinear coefficient d- kl {-2o),coM is often expressed in using 
Miller's delta [15] ^jkJ and the linear susceptibilities 

d ijk (-2 - e )Zji x„ Z kk « ijk 

= e o {n*(2o» - 1 ( {nfcaj ) • 1( {\(® ) • U^ lfc ■ <4) 

Here Xa^ 03 is the linear electric susceptibility at frequency 2co, and «j(2o3) is the index of 
refraction at 2o) both for light with the electric field polarized in the direction of the i-th 
principal axis of the crystal. The dispersion of <5^ has been demonstrated to be small, 
and for experimental purposes it is usually treated as a constant. 

Beam characterization measurements described in the next section show that the 
transverse distribution of our beam is nearly Gaussian in shape. Second harmonic 
conversion efficiency for cw beams with Gaussian transverse distribution in the low 
conversion limit is given by [11], [13] 

n =p /p = 2oj 2 d 2 „p l k (tc« 3 tr 0 c 3 ) , (5) 

'cw 2co io ett to u 

where p u and p 2(a are the powers of incident fundamental (after surface losses) and the 

output second harmonic (before surface losses), respectively. The Boyd and Kleinman 
focusing factor h(B,£) is a function of the walkoff parameter B and focusing parameter c, 

which are expressed by the formulas 

B = pJTkl 2 ( 6 > 

y and 

£ * / / (* w 0 2 ) • W 

Here p is the birefringent walkoff angle, w 0 is the 1/e amplitude radius at the beam waist, 
and k = n col c is the magnitude of the the fundamental w-avevector inside the crystal. 

For beam parameters used in these measurements, the approximations of lor 

the limiting case of weak focusing % « 1 are appropriate. With the added constraints ot 
optimum phasematching, no absorption, and birefringent walkoff aperture length less than 
crystal length / a = 7 t I/2 vv 0 /p < / , the focusing factor is given to the necessary accuracy by 
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h (B,4)»Z ( 1 - r/12 + / 4 Z 1 20 - r 6 / 1 344 + ••• ) - 

where t = 2B (2£ ) 1/2 . Since Eq. (8) assumes perfect phasematching, numerical solutions 
of the double integral defining h(B,t; ) [1 1] were used to calculate phasematching curves 
for angle tuned crystals. Equation (8) does not include bulk material losses, however 
using Eq (8) and approximating bulk losses as surface losses gave the required accuracy 
for a 1 cm KDP crystal. None of the other crystals used had significant absorption. 

The extension of equation (5) to pulsed harmonic conversion was accomplished with 
an integration over the observed pulse shape; u = j p(t ) dr = P At where u is the total 
energy, p(t ) is the instantaneous power, P is the peak power, and At is an effective 
pulse width. Here At w and At 2w w'ere obtained by numerical integration of the observed 
fundamental pulse shape and the square of the observed fundamental pulse shape 
respectively. This treatment is applicable when the pulse has no structure or modulation 
on the scale picoseconds or shorter which would make group velocity walkoff important. 
Rewriting equation (5) for pulsed harmonic conversion 


'calc 


Zco 


to 


P At 

2co 2co 

P At 
01 0) 


to . 2 J 3 

At it n £ n c 
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The subscript is used to indicate that 77 calc is the conversion efficiency that is obtained from 

a calculation that assumes no pump depletion. Pump depletion was approximated with 
the relation 

^calc ^observed ^ ^ ^ ^observed ^ ’ 


where ^observed ' s observed energy conversion efficiency. Numerical integration for 
the near field without walkoff showed equation ( 10) is accurate to better than 2% for 
^observed < 50% . 



Allowance was also made for an elliptical transverse distribution ot the pump beam, 
a condition that existed in some of these measurements. A pump traveling in the z- 
direction with intensity distribution which is described by 
L = ! o exp(-2r Vw^ 2 -2 y 2 /w y 2 ) 

is used. For the near field approximation, it is necessary to retain the walkoff dependance 
of the focusing parameter in the critical phase-matching direction which is chosen to be the 
r-direction. However, the beam area must be adjusted to allow for a different size in the 
noncritical y-direction. This is done with the substitution 

h (B4 ) -> (w x /w y )h (B4 x ) ; £ x = l/(kw^) . 


When the pump depletion approximation and the allowance for elliptical transverse 


distribution are included, equation (9) solved for becomes 
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( 12 ) 


In Eq. (12) both a>= 27 tc/A 0 and k = 2 Jt n / A 0 have been expressed in terms of the 
ffee-space wavelength Aq of the fundamental radiation. 

The expressions for the effective nonlinear coefficients of the materials used are 
reproduced in Table I. The effect of nonorthogonality of the extraordinary electnc field 
with the wavevector in the birefringent crystals is included. For negative uniaxial crystals 
this is done by replacing 0 with 0 + p ( 1 1 ) in expressions such as those given by 
Midwinter and Warner (16]. Here 9 is the phase-matching angle and p is the birefringent 
walkoff angle. For LilO^, which has a relatively large walkoff angle, the magnitude ot the 
value of obtained from a measured value of d t ff is reduced. Indices of retraction and 
phasematching parameters calculated from dispersion equations are given in Table II. 
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in. Experimental Setup 


Measurements were performed in sets that involved the comparison of two nonlinear 
crystals. One crystal was a well characterized reference, and the other was the crystal 
under test. Each set included absolute measurement of the individual crystals and a 
relative measurement. Phase-matching tuning curves were observed for both crystals, and 
surface reflections and total transmissions were measured at both fundamental and 
harmonic frequencies. Transverse beam characterizations were performed in the course of 
a measurement set, and temporal pulse shape was monitored throughout the 
measurements. Other calibrations and consistency checks were performed on a less 
regular basis. 

Figure 1 shows the main components of the experimental setup which consisted of a 
Q-switched Nd:YAG laser, spatial filter, and two beam arrangement in which the 
nonlinear crystals were measured and compared. An important feature of these 
measurements was the use of a single-mode injection-seeded laser for the 1.064 pm pump 
source. Temporal properties of the laser output were determined with a 0.4 ns rise time 
photodiode-oscilloscope combination and a Fabry-Perot interferometer with 0.03 cm - ' 
resolution. The measurements showed that the injection-seeded laser operated in a single 
axial mode. We expect that the spectral distribution was near the time- bandwidth limit for 
the 7.0 ± 0.2 ns full-width-at-half-maximum pulse. 

Pulse duration changed little over the course of these measurements. Typical 
oscillograms of both fundamental and harmonic pulsed are shown in Fig. 2. The leading 
and trailing edges of the fundamental pulse had exponential shapes with 1.4 ns time 
constant tor the rise and 3.0 ns for the fall. Numerical integration over the measured pulse 
shape yielded the effective pulse width At^ = 8. 1 ± 0.2 ns and the ratio At^JAi^ = 0.67 
±0.02. The ratio zl^co / = 8.3 x 10^ s'* ± 3% was used in evaluating Eq. (12). 

The time response and time base calibration of the photodiode-oscilloscope combination 
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were checked by observing the output of a mode-locked Nd:YAG laser. The time base 
was accurate to within 0.6%. 

Shot-to-shot energy fluctuations were caused by drifting in the alignment oi the 
spatial filter. Figure 3a shows a histogram of individual pulse energy distribution. Most 
of the harmonic generation measurements consisted of averages obtained with 100 pulses. 
Several such averages were used in a measurement. Using 100 shot averages simplified 
the problem of shot-to-shot energy fluctuations and was consistent with the knife-edge 
transverse beam characterizations which were done with several hundred pulses. The 
histogram of pulse energies shows reasonable agreement with a normal distribution tor 
which the probability of finding the fundamental energy between m w and k u + du^ is 
given by P(u J d where 

P(u ) =exp(-(w - u ) 2 / (2<7 2 )}/(<t - Jin ) , 

0) CO CO CO CO 

u is the average fundamental energy, and <7^ is the standard deviation. We observed 
a ju = 0.08. The standard deviation for harmonic pulses was 16%. However, the 100 

CO CO 

shot averages reduced this value by a lactor of ten. Averaging the shot-to-shot 
fluctuations will slightly bias harmonic measurements to higher values of harmonic power 
as is seen by calculating the average value of u ^ given by Eq. (13) 


(u 2 > = (u ) 2 + (7 2 

co CO CO 


The bias was only 0.6% and was not included in the calculation. 

Spatial filtering was performed by a combination of propagation into the far field, 
followed by aperturing the central peak ot the distribution, and finally focusing through a 
pin hole in a vacuum spatial filter followed by recollimation. The output beam was 
slightly converging and reached a beam waist about 2 meters beyond the collimating lens. 
Beam radii were measured by a knife edge method at several positions after the collimating 
lens. It was demonstrated that the beam propagation closely followed that predicted for a 
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Gaussian distribution. Aperture and pin hole adjustment in the spatial filter provided beam 
waists between w Q = 0.9 mm and 1.6 mm. 

Additional checks demonstrated that there were no extraneous temporal or spectral 
components in the laser output. No change in beam size due to a possible second spectral 
component was observed with the knife edge characterization when two prisms were used 
in place of turning mirrors to provide a dispersion of 4 nm/mm. Photodiodes with 
various integrating times and sensitivities were used with different triggering and sweep 
rate settings of the oscilloscope, and no observable temporal components other than the Q- 
switched pulse were found in the laser output. 

Repeated measurements of beam distribution were performed at the position where 
the SHG crystals were placed. The two beam arrangement allowed quick measurements 
of the 0.05, 0.1, 0.25, 0.5, 0.75, 0.9, and 0.95 transmission positions as the knife edge 
was scanned across the beam in one channel. The seven measurements were averaged to 
locate the center of the beam, and each of the six measurements other than the 0.5 position 
yielded a measurement of w x or w y . If the horizontal and vertical measurements agreed 
within experimental error they were combined to a single value w Q , and if not the beam 
would be treated as elliptical in shape. Figure 3b shows the beam transmittance for knife 
edge positions compared to values calculated for a Gaussian distribution Typical 
individual beam scans yielded ±3% standard deviation for waist measurements. Eight 
measurements made over the course of one day combined to yield a standard deviation of 
i5% which was taken as the accuracy to which the beam waist was known. 

The components ot the two-beam arrangement were easily rearranged and 
interchanged for a variety of measurements. The pump beam was transmitted through a 
thin film polarizer and a partially transmitting mirror to provide 2 - 4 mJ of linearly 
polarized 1.064 p.m radiation incident on the beamsplitter. Beamsplitting ratios from 
50/50 to 93/7 were used. The pyroelectric energy probes used to measure pulse energy 
had flat spectral response and could monitor both fundamental and harmonic. Absorbing 
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glass filters and highly reflecting dielectric mirrors were used to eliminate fundamental 
radiation and transmit harmonic radiation when necessary. Calibrated neutral density 
filters and diffusers were used to keep measured fundamental intensities within the range 
of the energy probes. In addition to absolute and relative harmonic generation 
measurements, the apparatus also was used to measure transmissions of crystals and 
filters. Surface reflectivities of the crystals at both fundamental and harmonic were 
measured with rearrangement of components, and transmission at 532 nm was measured 
with the harmonic output generated in a separate crystal. 

Three axes of rotational adjustment were used in the crystal mounts. Two of the 
axes were perpendicular to the direction of polarization and a third was parallel. A 
polarized alignment beam was used to determine whether the crystals w-ere properly 
oriented with respect to the pump beam. A half-wave plate was used to rotate the 
fundamental polarization for maximum harmonic generation for measurements with type- II 
crystals. The angular adjustments allowed rotation through a small range of azimuthal 
angle to check that the crystals were properly oriented for maximum All the crystals 
were properly oriented. However, it was necessary to use the 12 mm Bal^O^ crystal 
significantly off normal incidence to avoid problems of parallel-surface reflections. 

The calibration of the two pyroelectric energy probes were intercompared, and each 
pyroelectric probe was compared to a thermoelectric power meter. The power meter was 
internally calibrated with an electrical resistance heater. The 30 Hz repetition rate ot the 
laser allowed comparison of the average power measurement with the average pulse 
energy of the pyroelectric probes. The thermoelectric power meter read 3% higher than 
the higher of the two pyroelectric energy probes, and that probe gave reading 2 % higher 
than the other; all agreed within 5%. The relative difference of the two probes was 
retained in the analysis of data, and their average reading was used as the energy 
calibration. The reason for this choice was that conditions were more appropriate for the 
pyroelectric detectors, and the 5% difference is within the expected error of all the meters 
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and is too small to be significant. The 5% value was taken as the accuracy of the energy 
measurement, which also is the accuracy stated by the manufacturer. 

Energy measurement was also dependent on the accuracy of filter, attenuator, and 
beamsplitter calibration. Filters and attenuators were measured both in the two-beam setup 
described above and with a spectrophotometer. Comparison indicated accuracies of 1% in 
the measured transmissions. The beamsplitting ratio and probe sensitivity ratio were 
checked repeatedly to assure consistency in the measurements. Variations of 1/2 to 1% 
caused by differences in probe placement and pointing were found. Where possible 
critical parameters were measured with more than one method, and multiple or additional 
measurements were used for consistency checks and determination of accuracy. 

The overall accuracy ot these nonlinear coefficient measurements is estimated to be 
better than 10%. The average of the harmonic energy divided by the square of the 
fundamental energy was determined for each series of 100 shots in making 

absolute measurements. A further average was made of at least five and usually more 
such measurements of Typically the individual averages in one set of 

measurements would be consistent to ±2%. Larger variations were observed with 
repeated measurements that would entail recalibration, use of different samples of the same 
material, or different beam parameters. For example 12 measurements on three different 
samples of BaB 2 0 4 had a ±3.6% standard deviation for the value of d eff . Examination of 
Et}- (12), however, shows that can only be determined to the accuracy to which wq 
was know, and that was estimated to be 5%. The accuracy of quantities At lw / (At ) 2 
an d ** 20 )/ cou ^ double the uncertainty of the values of d^. 

IV. Measurements 

Twelve separate absolute measurements were made on three different BaB 2 0 4 
samples, more than used for the other materials. Also the other materials were all 
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measured relative to BaB 2 0 4 . This discussion starts with the barium metaborate 
measurements. The discussion is then expanded to include other materials both by 
absolute and relative measurement of the nonlinear coefficients. 

Two of the three BaB 2 0 4 samples were grown at the Stanford Center for Materials 
Research, and the third was grown at the Fujian Institute of Research on the Structure of 
Matter in the Peoples Republic of China. The Stanford Crystals were 4.1 and 1 1.9 mm 
long, and the Fujian crystal was 9.4 mm long. The tuning curves for the three BaB 2 0 4 
crystals agreed well with tuning curves obtained from dispersion equations and numerical 
evaluation of the double integral that defines the Boyd and Kleinman focusing factor 
h(B4). The calculated tuning curves were adjusted for pump depletion using Eq. (10). 
The observed and calculated curves for the 1 1.9 mm crystal are shown in Figs. 4a and b. 

The effective nonlinear coefficient obtained in these measurements was d eff (BaB 2 0 4 ) 
= 1.94±0.07 pm/V. The sign of the d 31 coefficient relative to that of d-> 2 is unknown tor 
BaB ? 0 4 , and no effort was made to determine the polarity of the crystals that were 
measured. Fortunately d 3 j is small compared to d 22 ; measurements of d 3 ^\ - 0.07 \dy^ 
[23] and W 31 l < 0.05 Idod [24] have been reported. Assuming that d 3] is negligible, the 
value ld 22 l = 2.16±0.08 pm/V is obtained from the measured d eff , calculated phase- 
matching parameters (Table I), and the appropriate definition of d c{{ (Table II). 

Expressing the result with two digits as ld 22 l = 2.2 pm/V more correctly conveys the 
accuracy of the measurement. 

Earlier measurements have placed ld 22 (BaB 2 0 4 )l = 4.1 d 3(3 (KDP) [23], [24], The 
result of 2.2 pm/V is 37% higher than these earlier results when d 36 (KDP) = 

0.39 pm/V is used. Even though ld 22 (BaB 2 0 4 )l = 5.7 d 36 (KDP) was observed here, 
these measurements nearly duplicated the 0.39 pm/V result tor KDP. 

Only a single KDP crystal was measured. The 10.4 mm long crystal was oriented 
for type-I phase matching, and surfaces were uncoated. The observed tuning curve again 
closely follows the calculated curve (Tigs. 4c and d). There is a smaller amount of 
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birefringent walkoff in this crystal compared to the BaB 2 0 4 crystal; therefore the tuning 
curves more closely approach the sine-squared curve of the monochromatic planewave. 
Both absolute measurements ot the KDP crystal and measurements relative to the 9.4 and 
1 1.7 mm BaB 2 0 4 crystals using d eff <BaB 2 0 4 ) = (1.9410.07) pm/V were closely grouped 
at d 36 (KDP) = 0 .37610.005 pm/V. This small range of measurements could be 
coincidence, and the earlier comments about accuracy apply. More typically the combined 
measurements on a single material would have ±4% standard deviation. For comparison it 
is useful to mention that in the notation used by Craxton [2] {</ 36 (KDP)/£^} Cra;aon = 0.78 
pm/V corresponds to ^-^(KDP) = 0.39 pm/V in the notation used here. 

In contrast to the agreement with KDP second harmonic measurements, these 
observations produced substantially different results compared to parametric fluorescence 
measurements of lithium iodate. Two LiI0 3 crystals 14.7 and 19.8 mm long were used 
tor type-I second harmonic generation, rhese crystals were antireflection coated on input 
and output faces. Both crystals demonstrated good agreement between predicted and 
observed phase-matching curves. Birefringent walkoff was again significant as can be 
seen in the shape of the secondary maxima in the tuning curves of the 14.7 mm crystal. 
(Figs. 4e and 0 Absolute measurements of the two LiI0 3 crystals gave lti 31 1 = 4.24±0. 10 
pm/V, whereas measurements relative to BaB 2 0 4 and KDP respectively yielded 
4.02±0.19 and 4.08±0. 19 pm/V. The value obtained by combining these measurements 
and using the known sign was </ 31 (LiI0 3 ) = -4. 1±0.2 pm/V. 

Tabulations of measurements of the nonlinear coefficients of LiI0 3 [6] show an 
unusually wide range of values. Of concern was that the high intensities used for pulsed 
harmonic generation could lower efficiency though some third order nonlinear effect such 
as two-photon absorption or intensity dependent refractive indices. Harmonic generation 
was observed as a function of fundamental pulse energy in an attempt to detect such 
effects. This was accomplished by attenuating the pump pulse with calibrated partially 
reflecting mirrors between the beamsplitter and the 19.8 mm LiI0 3 crystal. Figure 5a 
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shows the observed deviation from linearity between generated harmonic energy and 
square of the fundamental energy. The deviation form linearity was explained by pump 
pulse depletion as approximated by Eq. (10). There was no indication ot reduction of 
harmonic generation by any third order process in this set of measurements. The four 
steps of attenuation gave five relative measurements of nonlinear coefficient at peak 
intensities ranging from 7.4 MW/cm 2 for the unattenuated pulse to approximately one-fifth 
that value. The five measurements had better consistency, d 31 = -4.13±0.05 pm/V, than 
other measurements made at constant intensity. 

An additional measurement was performed to test for nonlinear absorption. 
Fundamental energy was monitored in the reference channel, and both transmitted 
fundamental and generated harmonic were measured in the test channel as the 19.8 mm 
LiI0 3 crystal was tuned through phasematching. The measurement of total transmission 
and observed and calculated tuning curves are shown in Figs. 6a, b, and c. The 
measurement suggested that two-photon absorption was present but not large enough to 
significantly have changed the measurement of nonlinear coefficient. The second surface 
of the crystal was antireflection coated for 532 nm, and a 0.5% increase in total 
transmission was expected at peak conversion but not observed suggesting a few percent 
of the harmonic was lost by two-photon absorption. 

Parametric fluorescence measurements have been previously performed tor lithium 
niobate [25], [7], and it was of interest to compare second harmonic generation 
measurements in this material. A 6.3 mm sample of 5%Mg0:LiNb0 3 was prepared for 
temperature-tuned noncritically-phase-matched second harmonic measurements. This 
crystal was uncoated. The previous parametric fluorescence measurements were 
performed with stoichiometric and congruent LiNb0 3 and not with the MgO doped 
material used for these measurements. The MgO doped LiNb0 3 has slightly different 
properties than those of the other compositions of LiNb0 3 but should be close enough tor 
a useful comparison of the two techniques ot nonlinear coefficient measurement. 
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The measured coefficient d 31 (5%Mg0:LiNb0 3 ) = - 4.7 pm/V is only 79% of the 
parametric florescence value for congruent LiNb0 3 , not unreasonable agreement 
considering that the materials have different composition. This is supported by 
comparison measurements [29] of second harmonic generation by 9 mm samples of MgO 
doped and congruent materials under similar pumping conditions; 50.9% conversion was 
observed in the congruent material, and 35.2% conversion was observed in the 
MgO:doped material. 

Two additional materials, KD*P and KTP, commonly used for 1064-nm-pumped 
harmonic generation were studied. Two 30 mm long KD*P crystals, one each of type-1 
and type-n were measured. For the type-II crystal, the result of a quantitative analysis 
derived from a near field calculation was used, and Kleinman symmetry was used to set 
d 36 = d \4- The KD*P crystals were measured both absolutely and relative to BaB-,0 4 
with the result d 36 (KD*P) = 0.367±0.012 pm/V. This is consistent with earlier 
measurements relative to KDP [30], [31], [32] which yielded values in the range 
d 36 (KD*P) = 0.34 to 0.48 pm/V when normalized to d 36 (KDP) = 0.39 pm/V [3] or 
^ 36 (KDP) = 0-41 pm/V [6], The tuning curves for the two KD*P crystal both indicated 
high quality crystals and little distortion of the pump beam. The tuning curve of the type-I 
crystal is shown in Figs. 4g and h. 

Second harmonic generation was observed in three type-II KTP crystals. All the 
KTP crystals were produced by the flux growth technique in the Peoples Republic of 
China. The very parallel surfaces of the first two crystals produced interference which 
was seen in the tuning curve (Fig. 7). The third crystal was also very parallel but had a 
1064 nm antireflection on one surface and a 532 nm antireflection coating on the second 
surface. The observed tuning curve for this crystal was in excellent agreement for 
rotations both about the z-axis and about an axis in the x-y plane as shown in Fig. 8. This 
agreement indicates a crystal of excellent optical quality, with proper orientation, and very 
small distortion of the pump beam. The polarization of the fundamental pulse was rotated 

- 19 - 


with a half-wave plate to produce the maximum second harmonic. The effective nonlinear 
optical coefficient derived from absolute harmonic measurements and measurements 
relative to KDP, KD*P, and BaB 2 0 4 was ti cff <KTP) = 3.18±0.17 pm/V. The value 
calculated from earlier measurements = 7.6xl0' 12 m/V and W 15 1 = 6.1 pm/V [33] is 

^KTPJcalc. = 7 ‘ 3 t 22] ‘ 

It is necessary to know the relative signs and relative magnitudes of d l5 and for 
the purpose of determining the size from a measurement of d cn . Observation of the 
change in d e ^ with crystal orientation indicates that d^ j and d -^ 2 have the same sign . If it 
is assumed that the ratio ^4^15 = 125 as measured previously [33], then the measured 
value of d^f yields — 3.3 pm/V and j = 2.6 pm/V. These values are less than half 
of the earlier values. The measurements, complicated by interference, made with the 
uncoated KTP crystals ranged from d cff = 2. 1 to 4.4 pm/V and supported the 
measurement of the lower value. 

V. Summary 

The nonlinear optical coefficient measurements reported here are summarized in 
Table III. These values can be taken as both relative measurements and absolute 
measurements. The reproducibility of the measurements was approximately ±4%. When 
the values are used as relative measurements, the accuracy is the same as the 
reproducibility. It is possible that the absolute measurements could be biased by 
inaccuracy in modeling and fluctuations of the transverse beam distribunon, measurement 
of pulse shape and duration, and by the accuracy of the energy measurements. It is 
estimated that the accuracy of the absolute measurements of the nonlinear coefficients is 
± 10 %. 

The measurement for KDP reported here is in agreement with the nonlinear 
coefficient accepted for that material in high power nonlinear conversion applications, and 
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it agrees with the combination of relative measurements between KDP and ADP and low 
power cw second harmonic measurements with ADP. Furthermore the relative 
measurements between KD*P and KDP and between KDP and LiI0 3 are in agreement 
with earlier measurements. There is possible agreement between the absolute 
measurement of 5%Mg0:LiNb0 3 and parametric fluorescence measurements of congruent 
LiNb0 3 , but this is not conclusive. We have measured a significantly lower value for the 
nonlinear coefficients of KTP and moderately higher value for the nonlinear coefficient of 
Bal^C^ than reported in earlier measurements. Perhaps the most significant difference is 
the absolute measurement of the nonlinear coefficient of LiI0 3 which is only 58% of the 
value measured by parametric fluorescence. 

It is remarkable, that after more than twenty-five years of study in nonlinear optics, 
that there should still exist such uncertainty in the scale of the nonlinear material 
parameters. This investigation indicates a need for further study. Performing both second 
harmonic and parametric fluorescence measurements on the same nonlinear crystal 
samples would provide useful information. The values for congruent for lithium niobate 
and MgO:LiNb0 3 will remain unresolved until this is done. The anomaly between second 
harmonic generation measured and parametric fluorescence measured nonlinear 
coefficients for lithium iodate also needs to be resolved. 

It is possible to make accurate measurements of nonlinear optical coefficients using 
the technique of second harmonic generation provided a great deal of care is used in 
making the measurements and high optical quality fundamental radiation is used. The 
development of automatic data acquisition and reduction techniques and the availability of 
well characterized highly coherent lasers will make further accurate measurements 
tractable. The wide range of nonlinear coefficient values that have been obtained over the 
past twenty-five years underscores the need to use well characterized high-quality sources 
of fundamental radiation in nonlinear optical frequency conversion. 
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Appendix 

Some detailed aspects of nonlinear frequency conversion in birefringent crystals are 
considered in this appendix. An anisotropic-media Green's function analysis is used to 
show that the treatment used here is appropriate for quantitative analysis of second 
harmonic generation under near field conditions. Only under tight focusing conditions will 
birefringence cause some astigmatism and aberration that is not already included in the 
analysis. 

A central assumption of the Boyd and Kleinman treatment of harmonic generation by 
focused beams [ 1 1] is that a Gaussian harmonic beam is generated in each infinitesimal slab 
increment of a nonlinear crystal through which a Gaussian fundamental beam is 
propagated. In the case of type-I phase matching in a negative uniaxial crystal, the 
fundamental beam is an ordinary wave and the harmonic is extraordinary. It is iurther 
assumed that the extraordinary harmonic beam has the same focal position / and the same 
confocal parameter b as the fundamental beam. However, the transverse position of the 
extraordinary beam is displaced in a way that causes it to propagate with the appropriate 
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birefringent walkoff. This result was derived by Kleinman, Ashkin and Boyd [34] using 
an isotropic-media Green s function modified to include birefringent walkoff. 

We proceed with the anisotropic-media analysis expressing results in a form similar to 
reference [34] to facilitate comparison. Assume that the nonlinear crystal exists between z 
= 0 and z = / t and is imbedded in a medium with identical linear birefringent properties but 
with no nonlinearity. The harmonic electric field at a point r is given by 

_ . exp(-2/ ox) I vr 

E 2 (r,t) = exp(-a (z-z')/2}YG(r,r>P‘ N V) dr' (A1 

4xe 0 J 

Here 2<y is the angular frequency of the harmonic, and (Xj is the absorption coefficient at 
the harmonic frequency. The Green s function G(r,r’) is described later. The dyadic 
projection operator is defined by 

4 Ki co A a 

Y = u u , (A2) 
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where n e (0) is the extraordinary index of refraction for k e (2to), and u is a unit vector in 
the direction of the electric field of the extraordinary wave specified by k e (2cu). The 
nonlinear polarization produced by an ordinary Gaussian fundamental beam propagating in 
the z- direction is 
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where x‘—2{z -f)/b and B (z ) — 1 tor G < z’ < 1, and 0 elsewhere. The fundamental 
wavevector is k x , and a, is the absorption coefficient of the fundamental. Everything to 
this point is identical with reference (34) and is presented for definition. 

The Green's function used here is correct for an anisotropic dielectric [10J 
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Equation (A4) differs from the Green's function used in reference |34] by the factor 
k K ^ cos p in the denominator, where K c is the Gaussian curvature ot the surtace 
defined by k e (2a>), and p is the birefringent walkoff angle. For an isotropic medium 

cos p is one, and the Green's functions become the same. As an example of an 
anisotropic material consider LilO^ oriented for phase-matched, type-I, 1064 — > 532 nm 
second harmonic generation tor which k I K cos p~ 1.10. The factor, however, cancels 


when the variation of k c with direction is included in the x' and y integration ot Eq. (A 1 ). 


The paraxial expansion used to perform the integration is 
exp{i k e (r— r ) } exp{2t k^(z - z )} 

I r — r' I (z - z' ) 
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Phase matching is assumed in this expression, and k x and K y are the principal curvatures 
of the surface defined by k c (note K c = k x K y ). Again for isotropic material we would have 
2&] = k^K x = 1 , and 2 k ] tC, = k Q K y = 1, and Eq. (A5) would become identical with the 

expansion used in ref. [34], 

Completing the x' and y' integrations of Eq. (Al) using (A4) and (A5) we have 

exp{2( k.z-a~z / 2) 

E (r) = r-P 0 x 

~ cos p 

f | -21 x-pjz-f_ )f_ _ 2y 2 J e~ az tfz' t (A6) 
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where T x = 2 { z - \f / k^K^- ( \/ k^K^ — \ )z ~\] k^K^/ b and (A7) 

T y = 2{z-\f/k e K y - ( l/k e K y - 1 )z’]}k c K y /b . (AH) 

We could funher manipulate the expressions (A7) and (A8) to have ^-distribution and y- 
distribution focal positions and confocal parameters, but it is unnecessary' for our purposes. 

We can project the solution back to the plane z - / , and take that as the exit surface ot the 
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crystal. With the experimental conditions of these measurements, we are in the near field 
with r« 1 at this plane. In using Eq. (8) we had assumed that x could be ignored, and 
the small changes in x given by the last two equations are of no significance. In the 
approximation that x x = r y = x= (z -f)/b, equation (A6) becomes the same as Eq. (4.19) 
of ref. [34], the result mentioned above which leads to Eq. (2.9) of the Boyd and Kleinman 
focused beam analysis [11]. 

For parametric fluorescence in negative uniaxial crystals such as LiI0 3 or LiNb0 3 , 
the scattered radiation or parametrically generated noise is ordinary, and an isotropic 
analysis of the scattering is adequate. The pump radiation used for these parametric 
fluorescence measurements, however, is extraordinary. Therefore it is necessary to use 
the angle of propagation plus the walkoff angle (9+ p) when calculating components of 
the extraordinary electric field of the pump radiation that lie in the direction of the principal 
axes of the crystal. In both parametric fluorescence measurements and second harmonic 
measurements in these crystals, the measured effective nonlinear coefficients must be 
converted to components of the nonlinear optical tensor using the angle (6 + p) [11], For 
LiNb0 3 with 9 = 90° , p =0, and there is no change from using 0 instead of (9 + p ). 
For LiI0 3 with 9 = 30.2° , p = 4.26° the effective nonlinear coefficient is increased by 
1.12 for a given value of d 3i using (9 + p ) instead of 9. This factor is not adequate to 
account for the discrepancy between measured here by SHG and previous 
measurements by parametric florescence. 
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Table I. Effective Nonlinear Optical Coefficients 


Crystal 

Point 

Group 

Phase- 

Matching 

Effective Nonlinear Optical Coefficient ^ 

KDP 

4 2m 

Type-I 

d eff = - d 36 sin (6 +p ) sin 20 

KD*P 

Type- II 

d e ff = (d l4 + d 36 ) sin (6+p ) cos (9+p) cos 20 

BaB 2 0 4 

5%Mg0:LiNb0 3 

3m 

Type-I 

d e[f = d 3l sin (6+p)- d 22 cos (6+p ) sin 30 

LiI0 3 

6 

Type-I 

d eff = d 3i sin i6+ P ) 

KTP0>) 

mm2 

Type-II 

d tf f = (d 24 - J 15 )sin20sin2p 
2 2 

- (d ] 5 sin 0+ d 24 cos 0 ) sin 6 


(a) There are different conventions used for defining effective nonlinear coefficients. We 
choose the one in which the direction of the wavevector is specified by sphencal 
coordinates (6,0 ) referenced to the crystalline axes; p is the birefringent walkotf 
angle. The positive sense of the extraordinary polarization is taken as that which has 
a component in the direction ot the positive crystalline z-axis.. 

(b) The expression for KTP is an approximation that is justified by the small difference 
between n T and n v compared to the difference between n z and n x or n , and further 
justified because 6 ~ 90° for all measurements used here. 
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Table II. Calculated Phase-matching Parameters 


Crystal Reference^ 

V 

P 

n 0 (co) 

n^co) 

n 0 (2co) 

« e (2a>) 

KDP, type-I 

[17] 

41.2° 

1.60° 

1.4942 

1.4603 

1.5129 

1.4709 

KD*P, type-I 

[18] 

36.6° 

1.45° 

1.4931 

1.4582 

1.5073 

1.4683 

KD*P, type-II 


53.7° 

1.42°, 1.28 

O 




BaB 2 0 4 , type-I 

[19] 

22.8° 

3.19° 

1.6545 

1.5392 

1.6742 

1.5547 

UIO3, type-I 

[20] 

30.2° 

4.26° 

1.8559 

1.7164 

1.8975 

1.7475 

5%Mg0:LiNb0 3 ^>, 

[21] 

90° 

0 

2.2327 

2.1527 

2.3242 

2.2327 

type-I 


T= 107 C 






KTP( C >, type-II 

[22] 

0 = 90° 

0.08®, 

n x ((0) = 

1.7381 

n x (2a)) = 

1.7785 



0 = 24.3° 

0.06° 

nJco) = 

1.7458 

n x (2co) = 

1.7892 





nlico) = 

1.8302 

n z (2co) = 

1.8894 


(a) Reference is to source of dispersion equations from which parameters are calculated. 


(b) LiNb0 3 is the only temperature tuned crystal; all the others are angle tuned at room 
temperature. 

(c) KTP is a biaxial crystal; all the others are uniaxial. 


4. 0 
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Table III. Nonlinear Optical Coefficients 


Crystal 

Nonlinear optical 
coefficient^ ( 1 O' 1 2 m/V) 

Miller's delta 
(10' 2 m 2 /C) 

KDP 

^36 = °' 38 

^36 = 2 - 4 

KD*P 

d 36 = 0.37 

5,6 = 2.4 

LiI0 3 

1 

II 

m 

■^3 

II 

1 

oo 

5%Mg0:LiNb0 3 

II 

1 

4^ 

5 31 = -0.84 

BaB 2 0 4 

d d{ - 1.94 
W 22 l = 2.2^ 

^t 

II 

"In 

KTP 

<W - 318 
W l5 1 = 2.6< c > 
1^24 1 35 3.3 

IS 15 I = 2.9 
l<5 24 l = 3.5 


(a) Nonlinear coefficients are given for 1064 to 532 nm second harmonic generation. 


(b) Assumes that ld 3 |l « ld 22 l for BaB 2 0 4 [24], 

(c) Using |d 24 l / |d 15 l = 1 .25 [33] and assuming d 2 4 and d {5 have the same sign. 
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Figure captions 


Fig. 1. Schematic drawing of experimental setup. 

Fig. 2. Oscillograms showing fundamental (a) and second harmonic (b) pulses. The 
time scale for both is 2 ns/cm displayed of the horizontal axis. Vertical displacement is 
relative power. 

Fig. 3. Measurements characterizing laser output, (a) Histogram shown fundamental 
pulse energy distribution, (b) Typical knife-edge measurement of transverse beam 
distribution and calculated transmission for a Gaussian distribution fit to the data 
points. Error bars were obtained using eight separate horizontal and vertical scans 
performed over the course of one day. 

Fig. 4. Measured (data points) and calculated (solid lines) phase-matching tuning curves 
for different second harmonic crystals. 

Fig. 5. Second harmonic generation energy conversion as a function of incident 
fundamental energy squared for lithium iodate (a) and lithium niobate (b) and (c). 
Fundamental energy was measured directly in the absolute measurement and was 
deduced for the harmonic generation in a reference crystal for the relative 
measurements. The solid lines are the calculated conversion with pump depletion, and 
the dotted lines without. Each data point represents a 100 shot average. The LiI0 3 
data points which fell obviously below the others were caused by unstable mounting 
of the crystal. 
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Fig. 6 . Total transmission for combined fundamental and harmonic as Li IO 3 (a) and 
LiNbOj (d) crystals are tuning through phasematching. The line drawn in the 
transmission measurement simply connects the data points to make them more visible. 
The measured phase-matching tuning curves are shown by data points and calculated 
tuning curves by the solid lines for LiI0 3 in (b) and (c) and for LiNb0 3 in (e) and (f). 

Fig. 7. A phase-matching tuning curve for an uncoated KTP crystal which has a problem 
of interference of reflections from the parallel surfaces. 

Fig. 8 . Phase-matching tuning curves of an antireflection coated KTP crystal tor 
noncruical rotation (a and b) and critical rotation (c and d). Data points show the 
observed tuning curves, and solid lines show the calculated tuning curves. 
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